The QCD phase diagram at nonzero baryon, isospin and strangeness chemical 
potentials: Results from a hadron resonance gas model 
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We use a hadron resonance gas model to study the QCD phase diagram at nonzero temperature, 
baryon, isospin and strangeness chemical potentials. We determine the temperature of the transition 
from the hadronic phase to the quark gluon plasma phase using two different methods. We find that 
the critical temperatures derived in both methods are in very good agreement. We find that the 
critical surface has a small curvature. We also find that the critical temperature's dependence on the 
baryon chemical potential at zero isospin chemical potential is almost identical to its dependence on 
the isospin chemical potential at vanishing baryon chemical potential. This result, which holds when 
the chemical potentials are small, supports recent lattice simulation studies. Finally, we find that 
at a given baryon chemical potential, the critical temperature is lowered as either the isospin or the 
strangeness chemical potential are increased. Therefore, in order to lower the critical temperature, 
it might be useful to use different isotopes in heavy ion collision experiments. 

Or 

C/5 ! I. INTRODUCTION 

i> : 

OA . The QCD phase diagram at nonzero temperature and baryon density has been the subject of numerous studies 
during the past few years 0,0- Two regions of the phase diagram have been placed under special scrutiny. First, in 
the low temperature and high density domain, the ground state is believed to be a color superconductor 
There are numerous types of color superconducting phases that emerge as the baryon chemical potential increases 
0,0. Second, in the high temperature and low baryon density domain, there is a transition from the hadronic phase 
^C) • to the quark gluon plasma phase Q. This latter domain is probed by heavy ion collision experiments. 

In heavy ion collision experiments baryon number, isospin and strangeness are conserved. The time between the 
formation of the fireball and its freeze-out is so short that only the strong interactions play a significant role. In the 
past heavy ion collision experiments, the value of the strangeness chemical potential is ~ 25% of the value of the 
baryon chemical potential, whereas the value of the isospin chemical potential is ~ 2% of the value of the baryon 
chemical potential (sj. Thus the strangeness and isospin chemical potentials are small but not negligible. It is 
Oh! therefore worthwhile to study the QCD phase diagram in the high temperature region with nonzero baryon, isospin 
ft, ■ and strangeness chemical potentials. 

| At zero chemical potentials, numerical simulations on the lattice find that the temperature of the transition from 
. the hadronic phase to the quark gluon plasma phase is given by T c = 175 ± 6 MeV [1 ©, Q3 El . However, at nonzero 
baryon chemical potential lattice simulations suffer from the so-called sign problem: The fermion determinant is 
• i-h , complex. As a consequence, traditional methods cannot be used to study the QCD phase diagram at nonzero baryon 
chemical potential. However, recent advances have allowed studies of the high temperature and low baryon chemical 
potential region 12, llMLahlOa. The situation at zero baryon and strangeness chemical potentials and nonzero 
isospin chemical potential is simpler: There is no sign problem and traditional methods can be used (as in QCD 
with two colors at nonzero baryon chemical potential jlSj - ). The results in this case are in complete agreement with 
effective theory studies 20] . From these lattice studies it appears that the critical temperature's dependence on small 
baryon chemical potential at zero isospin and strangeness chemical potentials matches its dependence on small isospin 
chemical potential at vanishing baryon and strangeness chemical potentials 0, fl3L ITiL Hq . We will show that the 
hadron resonance gas model supports this conclusion. 

Recently, the study of th e Q CD phase diagram at nonzero temperature, baryon and isospin chemical potentials 
has attracted attention pH l22l |23| . Using different models, it was found that an arbitrarily small isospin chemical 
potential could greatly alter the QCD phase diagram at high temperature and small baryon chemical potential and 
have important consequences for heavy ion collision experiments. There are two phase transitions at high temperature. 
There are phases where the u and d quark sectors are decoupled 0, |2^, [2j| . These results need to be confirmed by 
other methods, in particular on the lattice. 

We study the QCD phase diagram at high temperature and nonzero baryon, isospin and strangeness chemical 
potentials using the hadron resonance gas model. It has been shown both experimentally and on the lattice that the 
hadronic phase is very well described by a weakly interacting hadron resonance gas |gj, 124 , |25| . We use two different 
methods to determine the location of the transition. First, it has been found on the lattice that the phase transition 
that separates the hadronic phase from the quark gluon plasma phase corresponds to a surface of constant energy 
density: e ~ 0.5 — 1.0 GeV/fm 3 [25^. Second, the quark-antiquark condensate for the u and d quarks should almost 



2 



vanish at the transition temperature. In this article, we determine the critical temperature, T c , at nonzero baryon, 
isospin and strangeness chemical potentials by using both approaches. We compute the surfaces of constant energy 
density as well as the quark-antiquark condensate in a hadron resonance gas model at nonzero temperature, baryon, 
isospin and strangeness chemical potentials. We show that both methods agree qualitatively as well as quantitatively. 
We find that the critical surface has small curvature, and that the critical temperature slowly decreases when either 
the baryon, the isospin, or the strangeness chemical potentials are increased. 



II. THE MODEL 



We assume that the pressure in the hadronic phase is given by the contributions of all the hadron resonances up to 
2 GeV treated as a free gas, as in |2^. All the thermodynamic observables can be derived from the pressure since 

T 

P= lim 77 ln z ( t ^b,^i,^s,V), (1) 

V^oo V 

where Z(T, hb, [J-ii ^SiV) is the grand canonical partition function in a finite volume V, at nonzero temperature, T, 
baryon chemical potential, hb, isospin chemical potential, /i/, and strangeness chemical potential, us- The energy 
density is given by 

rj,dp dp dp dp 

oT d^ B <9m/ o^s 

For a quark q with mass m qj the quark-antiquark condensate is given by 

<«> = ^ q (3) 

At nonzero temperature, the contributions of massive states are exponentially suppressed ~ exp(— mn/T). Their 
interactions are also exponentially suppressed ~ exp(— (mn + m '-a)/T). Therefore this approximation should be valid 
at low enough temperatures. However, since we are studying QCD at temperatures up to ~ 200 MeV, the lightness of 
the pions could be a problem, since ~ 140 MeV. The hadron resonance gas model should be a good approximation 
for the other hadrons since they have a mass larger than ~ 500 MeV. The physics of pions at nonzero temperature 
has been extensively studied in chiral perturbation theory l26ll27j| . The pions' contributions to the pressure have been 
calculated up to three loops in chiral perturbation theory |26( . It has been shown that the free gas approximation and 
chiral perturbation theory agree at the one loop level, and that, in chiral perturbation theory, the two- loop corrections 
to the pressure are below a few percents of the one- loop contributions for temperatures under 200 MeV [2(|. Thus 
the hadron gas resonance model is a good approximation also for the pions. The hadron resonance gas model has 
already been used in the literature and has been shown to give a very good description of the hadronic phase and of 
the critical temperature |25j . 

In the free gas approximation, the contribution to the pressure due to a particle of mass mjj, baryon charge B, 
isospin .Z3, strangeness S, and degeneracy g is given by 

Ap = gm H T ^ -^-g- exp ^ j K 2 (— ) , (4) 

n— 1 v ' 

where 77 = +1 for fermions and r\ = — 1 for bosons, and K n [x) is the modified Bessel function. This particle's 
contribution to the energy density is given by 

^ - 9<T ± exp f ^s-I^-S^ ^ TR2 ^ + nmRKi ( « j (g) 

n=i n tt \ j 

and its contribution to the quark-antiquark condensate is given by 

A( qq ) = -gm H T — ^ exp j K, [—) . (6) 

In order to compute ©, we need to know tor as a function of m u or m^. We make two assumptions in order to 
compute ©. First, we assume that the Gell-Mann— Oakes— Renner relation is valid 

Fl ml = (m u + m d )(qq) Q , (7) 
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where F n — 93 MeV is the pion decay constant, and {qq)o — {uu)o — {dd)o is the quark-antiquark condensate at zero 
temperature and chemical potentials. Second, based on lattice results, we assume that the pion mass dependence of 
the hadron masses is given by 

d^H _ A 
d{mi) m H 

where 0.9 < A < 1.2 j^. Therefore, combining J2J and (JHJ), we assume that 

dmy ^ A(qq) Q 
dm q ~ F% mn ' 

Notice that since the hadron spectrum is isospin symmetric, we have that (uu) — {dd) = {qq) in the hadron resonance 
model at any temperature, bar yon, isospin and strangeness chemical potentials. Therefore, the rich structure of the 
phase diagram found in [2l], 122. |23| cannot be seen in this model. Finally at fixed T, it can be readily seen from (J5J 
and © that an increase in either hb, or fij, or /is will increase e and decrease {qq). Thus using either e or {qq) as a 
criterion to determine the critical temperature, we find that an increase in either fig, or fij, or fig decreases T c , and 
that at fixed /i^ an increase in \i\ or (is results in a decrease of T c as well. 



III. RESULTS 



A. Energy density criterion 

Lattice simulations have shown that the transition from the hadronic phase to the quark gluon plasma phase takes 
place at a constant energy density e ~ 0.5 — 1.0 GeV/fm 3 |25j. We use this criterion to determine the critical temper- 
ature as a function of baryon, isospin and strangeness chemical potentials. Our results for the critical temperature 
as a function of \jlb at fixed \ij and /is are shown in Figure ^ We find that this criterion constrains the critical 
temperature in a band of ~ 15 MeV. At zero chemical potentials, we find that T c = 176 ± 8 MeV, which is in good 
agreement with lattice simulations The temperature decreases as hb increases, as expected, but the 

decrease is slow. At the accuracy we can achieve using this method, an increase in /ij does indeed decrease the critical 
temperature at fixed /z#, but this effect is small at best. The decrease of the critical temperature is more important 
when fig is increased. In Figure [21 we compare the critical temperature as a function of (is at A*/ = Ms = with the 
critical temperature as a function of /i/ at (jlb = = 0. Notice that we limit ourselves to /i/ < and /is ^5 ffiK 
in order to avoid the pion and kaon superfluid phases 0, |2(], 13 ■ We find that the critical temperature curves are 
almost identical in both cases. This is in agreement with results from the lattice [T^. Il3l 113. Il5l Il8l| . 

We can fit our result for the critical temperature as a function of /is, \xi and /is- By construction, since the pressure 
is an even function of /i/ in the hadron resonance model, the critical temperature is also even in /x/. We find 

^ = 1 - 0.021(2) {^j 2 - 0.039(1) (^) 2 - 0.037(2) f^j - 0.031(3) + • ■ • , (10) 

where To is the critical temperature at zero chemical potentials. The fit is excellent, with a linear regression coefficient 
R 2 = 0.994. 



B. Quark-antiquark condensate criterion 

The critical temperature can also be computed from the quark-antiquark condensate. Indeed {qq) is of the order of 
the light quark masses at the phase transition and therefore almost vanishes. We determine the critical temperature 
by finding the point where {qq) — in the hadron resonance gas model. We obtain a range of critical temperatures, 
since in the relation 10 the constant A ~ 0.9 — 1.2 is not precisely known 25]. Our results for the critical temperature 
as a function of \xb at fixed fii and /is using {qq) are shown in Figure We find that this criterion constrains the 
critical temperature to a band of ~ 15 MeV. At zero chemical potentials, we find that T c = 185 ± 6 MeV, which 
is in good agreement both with the result obtained above using the e-criterion, as well as with lattice simulations 
H, 0, [nj El As expected, the critical temperature decreases as /is increases. As in the previous method, we find that 
an increase in the isospin chemical potential might reduce the critical temperature, but not in a significant way. The 
strangeness chemical potential has a stronger effect on the critical temperature. In Figure we compare the critical 
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FIG. 1: Critical temperature as a function of /is determined by lines of constant energy density: e ~ 0.5 — 1.0 GeV/fm 3 . In 
the upper two plots /is = 0, the dark shading with full curves corresponds to fii — 0, and the light shading with dashed curves 
corresponds to fii — 100 MeV. In the lower two plots [ii = 0, the dark shading with full curves corresponds to /is = 0, and the 
light shading with dashed curves corresponds to us = 200 MeV. To is the critical temperature at zero chemical potentials. 
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FIG. 2: Critical temperature as a function of /is at /i/ = /is = (dark shading with full curves), and as a function of /i/ at 
fiB = l^s = (light shading with dashed curves), determined by lines of constant energy density: e ~ 0.5 — 1.0 GeV/fm 3 . To is 
the critical temperature at zero chemical potentials. 



temperature as a function of \ib at \ii = /ig = with the critical temperature as a function of \ij at \ib = Ms = 0. 
We find that the critical temperature curves are almost identical in both cases. 

We can fit our result for the critical temperature as a function of /j,b, Hi and fis- We find 

Y q = 1 - 0.017(1) (^j 2 - 0.109(4) (£0 2 - 0.032(2) (^j ' - 0.024(2) + • • • , (11) 

where To is the critical temperature at zero chemical potentials. The fit is excellent, with a linear regression coefficient 
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FIG. 3: Critical temperature as a function of n q determined by (qq) = 0. In the upper two plots us = 0, the dark shading with 
full curves corresponds to hi — 0, and the light shading with dashed curves corresponds to fii = 100 MeV. In the lower two 
plots fii = 0, the dark shading with full curves corresponds to (is = 0, and the light shading with dashed curves corresponds 
to us = 200 MeV. To is the critical temperature at zero chemical potentials. 



T (MeV) 



(qq) 



Hb {hi = A*s = 0) 



Hi (mb = Ms = 0) 



T/T !- c 



0.96 



(qq) 



Hb (mj = Ms = 0) 



V 

Hi (Hb = Hs = 0) \ 



100 

H (MeV) 



0.2 



0.4 



FIG. 4: Critical temperature as a function of (xb at (ii = (is 
at hb — Hs — (light shading with dashed curves), determined by (qq) 
potentials. 
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R 2 = 0.991. 

Finally we can compare the critical temperatures obtained using these two different approaches: constant energy 
density and disappearance of the quark-antiquark condensate. We present our results in Figure[S]at nonzero hb with 
fij = /is = 0, nonzero /!/ with hb = Hs = 0, and nonzero mus with ^tg = \ij = 0, respectively. We find that 
the critical temperature at zero chemical potentials, To, is lower when we use the e-criterion than when we use the 
(gq)-criterion. If the critical curve is normalized with respect to To, we find that the two methods are in very good 
agreement. If we compare the fits ijTU)) and (fTTfl . we find that the hb and hs coefficients are very close in both cases, 
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FIG. 5: Comparison of the critical temperatures as a function of \xb at fii = and /is = (upper panel), as a function of /x/ 
at = and us — (lower left panel), and as a function of /ig at [is = and /i/ = (lower right panel) obtained using the 
energy density method (dark shading with full curves) and the quark-antiquark condensate method (light shading with dashed 
curves). Tq is the critical temperature at zero chemical potentials. 



whereas the \ij coefficients almost differ by a factor of three. However, this large difference in the coefficients leads 
to critical temperatures that only differ by a few percents in the region of interest. We therefore conclude that both 
methods yield critical temperatures that are in very good agreement. 



IV. CONCLUSION 



We have used the hadron resonance gas model to determine the temperature of the transition from the hadronic 
phase to the quark gluon plasma phase as a function of baryon, isospin and strangeness chemical potentials. This is 
of interest for heavy ion collision experiments, since baryon, isospin and strangeness chemical potentials are nonzero 
in this case. We have used two different methods to determine the critical temperature. The first one relies on the 
observation on the lattice that the quark gluon plasma phase emerges at a constant energy density [25|. The second 
method is based on the fact that the quark-antiquark condensate for the light quarks should almost disappear at the 
transition between the hadronic phase and the quark gluon plasma phase. We find that the critical temperatures 
found in both methods are in very good agreement. 

In the hadron resonance gas model, the critical temperature decreases as the baryon, isospin, or strangeness chemical 
potentials increase, albeit slowly. In agreement with recent lattice simulations [El Ha . X\M Hoi Il8 | and several models 
we find that the critical temperature as a function of the quark chemical potential at zero isospin chemical 
potential is almost identical to the critical temperature as a function of isospin chemical potential at zero quark 
chemical potential. We also find that the critical temperature decreases slightly when the isospin chemical potential is 
increased at fixed baryon chemical potential. This might be important for heavy ion collision experiments: A choice 
of different isotopes should reduce the critical temperature that separates the hadronic phase from the quark gluon 
plasma phase. 
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